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ABSTRACT 

The expulsion of axisymmetric magnetic field from the event horizons of rapidly ro- 
tating black holes has been seen as an astrophysically important effect that may sig- 
nificantly reduce or even nullify the efficiency of the Blandford-Znajek mechanism 
of powering the relativistic jets in Active Galactic Nuclei and Gamma Ray Bursts. 
However, this Meissner-like effect is seen in vacuum solutions of black hole electro- 
dynamics whereas the Blandford-Znajek mechanism is concerned with plasma-filled 
magnetospheres. In this paper we argue that conductivity dramatically changes the 
properties of axisymmetric electromagnetic solutions - even for a maximally rotat- 
ing Kerr black hole the magnetic field is pulled inside the event horizon. Moreover, 
the conditions resulting in outgoing Poynting ffux in the Blandford-Znajek mecha- 
nism exist not on the event horizon but everywhere within the black hole ergosphere. 
Thus, the Meissner effect is unlikely to be of interest in astrophysics of black holes, at 
least not in the way this has been suggested so far. These conclusions are supported 
by the results of time-dependent numerical simulations with three different computer 
codes. The test problems involve black holes with the rotation parameter ranging 
from a = 0.999 to a = 1. The pure electrodynamic simulations deal with the structure 
of conductive magnetospheres of black holes placed in a uniform-at-infinity magnetic 
field (Wald's problem) and the magnetohydrodynamic simulations are used to study 
the magnetospheres arising in the problem of disk accretion. 
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1 INTRODUCTION 

The electromagnetic mechanism of extraction of rotational 
energy of black holes by Blandford and Znajek l|1977p is con- 
sidered as one of the most promising models for powering 
the relativistic jets of active galactic nuclei (AGN) , gamma- 
ray bursts (GRBs), and galactic black hole binaries. When 
this mechanism is described in terms of classical physics the 
black hole horizon is often compared with a rotating mag- 
netized conductor (|Damour 19781 IRuffini fc Treves 1973p . 
This viewpoint was canonized in the "Membrane paradigm" 
by Thorne et al. I|1986p . Such a description puts black holes 
on the same footing as say the Sun since the Blandford- 
Znajek mechanism appears not much different from the 
mechanism of magnetic braking by Weber & Davis (|1967|> . In 
this respect the so-called "Meissner effect of black hole elec- 
trodynamics," that is the expulsion of magnetic ffux from 
the horizon of rapidly rotating black holes, seems to under- 
mine the role of the Blandford-Znajek mechanism in astro- 
physics as "the conductor" becomes unmagnetized. 

Here is the brief history of the Meissner effect. 



Wald l|1974p found the exact steady-state vacuum solution 
for a rotating black hole placed into a uniform-at-infinity 
magnetic field aligned with hole's rotational axis. For this 
solution. King et al. (|1975| ) computed the total magnetic fiux, 
$, threading the event horizon as a function of hole's rota- 
tion parameter, ^ a < 1. They found that <E> ^ as 
a 1 and, thus, for maximally rotating black holes, the 
magnetic fiux is totally expelled from the black hole horizon 
(see fig. By its appearance the phenomenon is similar to 
the Meissner effect, that is the expulsion of magnetic field 
by superconductors. Later Bicak & Janis (|1985p showed that 
this result held for all axisymmetric steady-state vacuum so- 
lutions and concluded that their finding has a detrimental 
effect on the prospects of the Blandford-Znajek mechanism. 
Last year this argument was reiterated by Bicak et al. (|2006p . 

Recently, there has been impressive progress in 
numerical methods for general relativistic magneto- 
hydrodynamics l|Koide et al. 19991 IKomissarov 20011 

De Villiers fc Hawley 2003] IGamm ie et al. 20031 

Komissarov 2004bl IDuez et al. 20051 ^nton ot al. 2006] 
Anninos et al. 20061 [Shibata fc Sekuguchi 2005] 
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lAnderson et al. 2006|l and a number of groups have 
carried out simulations of magnetized flows around Kerr 
black holes. None has reported observations of the Meissner- 
like effect and so it is natural to ask why. McKinney & 
Gammie l|2004p suggested that this could be due to a 
relatively small value of parameter a in their simulations - 
indeed, the effect is not strong unless a is very close to unity. 
This explanation has also been put forward in Bicak et 
al.(. 2006J . However there are other reasons that may be even 
more important. Indeed, the difference between the systems 
of magnetohydrodynamic (MHD) and vacuum electrody- 
namics is very significant both in structure and number 
of equations. The two physical factors that are accounted 
for in MHD but not in vacuum electrodynamics are the 
conductivity and inertia of plasma. Both of them have been 
argued to be of crucial importance for powering relativistic 
outflows from black holes ( jSlandford fc Znajek 1977] 
[Punsly fc Coroniti 1990"| ). 

In this paper we show that conductivity alone is suffi- 
cient to nullify the Meissner effect even for maximally rotat- 
ing black holes. In particular, we present the results of elec- 
trodynamic simulations for the Kerr black hole with the ro- 
tation parameter a = f .0 placed in an aligned uniform mag- 
netic field. We also present the results of a general relativistic 
magnetohydrodynamic (GRMHD) simulation of a realistic 
magnetized accretion flow with black hole spin a = 0.999. 
The GRMHD model confirmes that the inclusion of iner- 
tia and thermal effects does not help to revive the Meissner 
effect. 



2 ELECTRODYNAMICS OF BLACK HOLES 

In this study we only consider the astrophysically most rele- 
vant case of test fields, that is we assume that the mass- 
energy of the electromagnetic field is too small to effect 
the curvature of space-time. In our analysis we employ the 
3-f 1 formulation of black hole electrodynamics developed in 
Komissarov (|2004ap . While the better known formulation 
of Macdonald and Thorne (1982) is adapted to the Boyer- 
Lindquist foliation of spacetime that introduces a coordinate 
singularity on the event horizon, the formulation in Komis- 
sarov (|2004aP is more general and can be used for the Kerr- 
Schild foliation as well. This is important because there is 
no coordinate singularity in the Kerr-Schild foliation, which 
therefore has advantages for both numerical and analytical 
studies of physical processes in the vicinity of the event hori- 
zon. In particular, this allows us to place the inner boundary 
of the computational domain inside the event horizon. The 
differential equations of this formulation are identical to the 
equation of classical electrodynamics in matter 



V B = 0, 

9tB + V X £ = 0, 

V D = p, 

dtD + VxH ^ J. 



(1) 
(2) 

(3) 

(4) 



The purely spacial vectors B,D,E,H axe related to the 
Maxwell tensor via 



2 ^ 



(5) 
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where Cijk is the Levi-Civita tensor of space. 

The curved space of black hole behaves as an electro- 
magnetically active medium whose electromagnetic proper- 
ties are described by the following constitutive equations 



E ^ aD + fBxB, 
H = aB - l3xD, 



(9) 
(10) 



where a is the lapse function and /3 is the shift vector of 
the spacetime foliation. Following the general principle of 
relativity it is most convenient to describe local physical 
processes in the frames of local inertial observers who are 
instantaneously at rest in the "absolute space" of the folia- 
tion ("fiducial observers" or FIDOs). Such an observer sees 
B as the local magnetic field, D as the local electric field, 
and p as the local electric charge density. The local electric 
current density, j, is related to J of the Ampere equation 



J ^ aj - pf3. 



(11) 



Neglecting the inertia of charged particle one can write 
the generalized Ohm law as 



where 



3d 



DxB 
J — = 

£2 



(12) 



(13) 



is the drift current, is the conductivity current parallel 

to the magnetic field and a±D± is the conductivity current 
perpendicular to the magnetic field l|Komissarov 2004ap . In 
a coUisionless plasma with a strong magnetic field the cross- 
field conductivity is highly suppressed so one can safely use 
a± = 0. On the contrary, the parallel conductivity is very 
large, thus resulting in almost vanishing parallel component 
of the electric field, Dm <^ B. 



3 COMPUTER SIMULATIONS 

In our search for insights into the role of conductivity in 
the Meissner effect of black hole electrodynamics, we first 
consider the famous Wald problem (|Wald 1974[> . The steady 
state solution of vacuum electrodynamics to this problem 
describes not only the magnetic field but the electric field as 
well. The magnetic part is 



B"" = ~Bo{g^4,,B + 2agt,p,B)/y/^, 
B" = Boi 



(14) 



+ 2agt^^r)/y/l, 



where t, 0, r, 6 are the Kerr-Schild coordinates, g^^ is the 
metric tensor of spacetime, and 7 is the determinant of the 
metric tensor of space. The electric part is 



E, 



= 0, 



Er — 
Eg = 



-Bo{ 
-Boi 



-jt4,,r + 2ag 
htl>,e + 2agtt,0 



r), 



(15) 



'Meissner effect" and Blandford-Znajek mechanism in conductive black hole magneto spheres 3 




-2 2 -2 2 



X X 

Figure 1. Left panel: Steady-state numerical solution for the vacuum case. The magnetic flux is expelled from the event horizon 
thus illustrating the "Meissner effect" of black hole electrodynamics. Right panel: Numerical solution for the conductive case at time 
t = dGM/c? . The magnetic flux is pulled back onto the horizon. The unit length corresponds to GM/c^ . The thick line shows the black 
hole ergosphere (a = 1.0). 




Figure 2. Left panel: Similar to the right panel of figure [T] showing the steady-state numerical solution for the highly conducting case. 
All of the flux passing through the ergosphere threads the black hole. Right panel: Same model but shows field angular frequency of 
rotation, where for x > there are 11 contours from Q,p = to Vlp = Vljj/2. This demonstrates presence of cylindrically coUimated 
Poynting jet. 



(The components of B and E are given in the non- 
normalized coordinate bases of Kerr-Schild coordinates.) 
The emergence of the electric component in this solution is a 
particular example of the electromagnetic activity of curved 
space-times. This effect is of great astrophysical importance 
as it leads to the Blandford-Znajek process. 

The equations of conductive electrodynamics of black 
holes are more complicated than those of vacuum electro- 



dynamics and the analytic approach has not been very suc- 
cessful so far. Fortunately, there are now a number of nu- 
merical techniques that one may use. In our study we first 
used the upwind scheme described in Komissarov (2004a)) 
and set (T|| = 1/At where At is the computational time- 
step. Since At is much smaller than the hght crossing time 
of the magnetosphere, the solution is expected to be fairly 
close to one with infinite conductivity, and this is confirmed 
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by the results. We used solution (|14ll5p as the initial so- 
lution of our numerical simulations and to make sure that 
the Meissner effect is apparent we set a = 1.0, that corre- 
sponds to a maximally rotating black hole. We carried out 
two runs: one with vanishing conductivity, cr|| = 0, and the 
other with cr|| = 1/ At, which corresponds to the highly con- 
ductive case (except in current sheets). The left panel of 
figure [1] shows the steady-state numerical solution of vac- 
uum equations. One can see that magnetic flux is expelled 
from the black hole horizon. In fact, on this plot the nu- 
merical solution is indistinguishable from the exact analytic 
solution. The right panel of figure [T] shows the conductive 
case at time t — 6GM/c^. Now the magnetic flux is no 
longer expelled from the event horizon but on the contrary 
the magnetic fleld lines are actually attracted to it. By this 
time the solution starts to develope the dissipative equa- 
torial current sheet (|Komissarov 200 4a|. Unfortunately, the 
code cannot handle this kind of current sheets where the 
electric field tends to become stronger than the magnetic 
field. It turns out that the numerical structure of such cur- 
rent sheets in this code is controlled by numerical resistivity 
HKomissarov 2006il . Thus this scheme gives reliable results 
only until the ergospheric current sheet is formed. 

McKinney (|2U06a[) constructed a different numerical 
scheme that provides greater control over the influx of en- 
ergy into such currents sheets thus allowing them to evolve 
to an almosts dissipationless state. Figure [5] shows the al- 
most steady-state solution to the Wald problem with a = 1.0 
obtained with this code (a low level unsteady reconnection 
controlled by numerical resistivity occurs in the equatorial 
current sheet.) The left panel of flgure[2]shows the magnetic 
field lines and the right panel shows their angular velocity , 
Q,F = —Eg/y/^B'^, where 7 is the determinant of the met- 
ric tensor of space. This result is similar to that obtained 
by Komissarov (f2005)) for the GRMHD Wald problem with 
a = 0.9 where the reconnection rate in the equatorial cur- 
rent sheet is inhibited by the high thermal pressure of locally 
heated plasma. The Blandford-Znajek power from the black 
hole is maximally efficient since much of the magnetosphere 
rotates with the angular velocity S7f = Q.n/'i, where 
is the angular velocity of the blach hole. The cylindrically 
collimated Poynting jet is driven by both the black hole and 
the ergospheric disc. 

Finally, we also use the GRMHD code HARM 
UGammie et al. 20031 [McKinney 2006b| to confirm the re- 
sults of electrodynamic simulations, to see if the inclusion of 
plasma inertia and thermal pressure makes any difference, 
and to study the problem in the astrophysically more real- 
istic context of disc accretion. Here we examine a GRMHD 
solution of a magnetized accretion onto a black hole with 
a = 0.999 and otherwise similar parameters as described in 
McKinney & Gammie (pOOi)) . Figure [3] shows the GRMHD 
solution for the field geometry and rest-mass density. We 
find that, as in the previous simulations with a slower ro- 
tating black hole ( [McKinney 2005] ), there is no sign of mag- 
netic fiux expulsion from the event horizon. In the funnel 
the magnetic geometry is close to the monopole assumed 
in the model by Blandford and Znajek (1977). Thus, even 
for very rapidly rotating black holes the Blandford-Znajek 
effect can drive magnetized jets similar to those seen in sim- 
ulations with slower rotating black holes ( [McKinney 2006bl 
[Hawley fc Krolik 2006[ ). 




. Contours show the magnetic field lines and the colour image 
shows logarithm of the rest-mass density. Amount of flux thread- 
ing black hole is qualitatively similar to lower black hole spins for 
similar models. The Meissner efi'ect is not observed. 



4 DISCUSSION AND CONCLUSIONS 

The numerical simulations described here show that the 
"Meissner effect" is nullified when high conductivity, typ- 
ical for the magnetospheres of astrophysical black holes, is 
taken into account. What are the reasons for this result and 
how general is it? In order to answer these questions it is 
perhaps helpful to speculate on why one would expect the 
results of vacuum electrodynamics to be applicable under 
the conditions where conductivity is undoubtedly very high. 
Consider for example the case of a uniform magnetic field 
in flat spacetime. This is an exact steady-state solution of 
vacuum Maxwell's equations. Now lets us uniformly fill the 
space with static electrically neutral plasma. Obviously, this 
does not disturb the equilibrium and the magnetic configura- 
tion remains unchanged. The same applies to any other equi- 
librium magnetostatic configuration (we ignore gravity.). 

The difference between this case and the case of ro- 
tating black holes is in the fact that strong electric field is 
necessarily induced in the vicinity of a black hole when it is 
placed in vacuum magnetic field even if the black hole itself 
has zero electric charge. In the particular case of the Wald 
problem this "gravito-rotationally induced" electric field is 
described by equations (|9ll5p . However, this result is quite 
generic ([Komissarov 200"4a|) . Once plasma is introduced in 
such an electromagnetic field the charged particles of differ- 
ent signs will move in the opposite directions bringing about 
screening of the electric component. This screening deeply 
upsets the equilibrium of vacuum solution as the electric and 
magnetic fields are tightly coupled via equations (|9I10|) . The 
effect is strongest in the vicinity of the event horizon where 
the magnitudes of electric and magnetic fields are compa- 
rable. Moreover, it has been shown in Komissarov (|2004a|) 
that no static distribution of electric charge can give full 
screening of the electric component within the black hole 
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ergosphere where under this condition the purely poloidal 
magnetic field is bound to have lower magnitude than the 
normal component of the electric field. To achieve marginal 
screening a poloidal electric current must flow through the 
ergospheric region thus strengthening the magnetic field by 
creating the azimuthal component that was absent in the 
vacuum solution. The hoop stress associated with this com- 
ponent tends to pull the magnetic flux back on to the event 
horizon. 

The disscussion of the Meissner effect with connection 
to the Blandford-Znajek mechanism in Bicak et al. (|2006p 
highlights the widely spread misinterpretation of the electro- 
magnetic mechanism by Blandford-Znajek, namely that the 
energy is generated by the rotating event horizon. In fact, it 
is not the rotation of the horizon, or the "stretched horizon" 
HThorne et al. 1986p . that makes possible the electromag- 
netic extraction of the rotational energy of black holes. It is 
much more instructive to consider the space itself as rotat- 
ing and this rotation is particularly strong within the black 
hole ergosphere where it results in unavoidable rotation of 
plasma in the same sense as the space (or the black hole). 
Punsly and Coroniti ( 1990 bl were the first who correctly 
and clearly indentifled the ergosphere as the key region for 
the magnetic energy extraction and constructed a plausable 
model of "inertially driven wind" based on this plasma ro- 
tation within the ergosphere. However, in their critisim of 
the Blandford-Znajek mechanism they failed to see that this 
mechanism is also based on the extreme conditions existing 
within the ergosphere and not on the properties of the event 
horizon. The signature of space rotation, expressed by the 
shift vector /3, is present outside of the event horizon, where 
it provides the additional coupling of the electric and mag- 
netic fields that holds even in steady-state configurations 
(see eas l9ll0| l. Screening of the gravito-rotationally induced 
electric field, that is the field generated via this additional 
coupling, involves generation of stationary poloidal electric 
current (|Komissarov 2004ap and this current remains non- 
vanishing even when a steady-state is reached. It is this cur- 
rent, driven by the marginally screened ergospheric electric 
field, that provides for the energy and angular momentum 
extraction in the Blandford-Znajek mechanism and as far as 
the magnetic fiux is not expelled from the ergosphere, and it 
is not expelled even in the vacuum solutions, the electromag- 
netic mechanism will continue to operate. Thus, even if the 
Meissner effect did hold in the conductive regime this would 
not be very important for the Blandford-Znajek mechanism. 
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